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Abstract 

We show that in the case of smah R-parity and lepton number breaking couphngs, 
primordial nucleosynthesis, thermal leptogenesis and gravitino dark matter are 
naturally consistent for gravitino masses ^ 5 GeV. We present a model 

where R-parity breaking is tied to B-L breaking, which predicts the needed small 
couplings. The metastable next-to-lightest superparticle has a decay length that is 
typically larger than a few centimeters, with characteristic signatures at the LHC. 
The photon flux produced by relic gravitino decays may be part of the apparent 
excess in the extragalactic diffuse gamma-ray flux obtained from the EGRET data 
for a gravitino mass m3/2 ~ 10 GeV. In this case, a clear signal can be expected 
from GLAST in the near future. 



1 Introduction 



Most supersymmetric extensions of the standard model impose R-parity [T] as an exact 
symmetry of the supergravity Lagrangian. In this way, one forbids renormahzable baryon 
and lepton number violating interactions which might cause too rapid proton decay [2]. 
On theoretical grounds, however, theories with and without R-parity are on the same 
footing, and in low-energy effective theories obtained from string compactifications R- 
parity plays no preferred role. 

One can also construct supersymmetric extensions of the standard model without 
R-parity [3], and the phenomenological constraints on these theories have been studied 
in great detail P]. Without R-parity conservation, the lightest superparticle (LSP) is no 
longer stable and, in general, it does not contribute to dark matter. 

Stringent constraints on the lepton number and R-parity violating interactions 

W^L=i = Kkjke'jlk + KjidiQjlk (1) 

are imposed by baryogenesis. Both operators contain lepton doublets. Together with 
sphaleron processes they therefore influence the baryon asymmetry at high temperature 
in the early universe. The requirement that an existing baryon asjTumetry is not erased 
before the electroweak transition typically implies [5] 

A,A'<10~^ (2) 

It is very remarkable that for such a small breaking of R-parity a gravitino LSP has a 
lifetime much longer than the age of the universe [6] . This is due to the double suppression 
by the inverse Planck mass and the R-parity breaking coupling, oc A^m^yg/^p^- We 
find for the gravitino lifetime 

.,,.,0^.. (_^)-', ,3) 

which is consistent with gravitino dark matter. 

For a gravitino LSP, the properties of the next-to-lightest superparticle (NLSP) are 
strongly constrained by primordial nucleosynthesis (BBN). In the particularly interesting 
case of a charged NLSP, like a scalar r- lepton, its lifetime has to be relatively short, 
''"NLSP ^ 10^ — 10^ s [7p, which typically requires 1713/2 < 1 GeV. Even for neutral 
particles, BBN excludes a neutralino NLSP for lifetimes longer than 10^ s due to the 



^See also [8]. Here, we consider mNLSP = 0(100 GeV). For a heavier charged NLSP, tonlsp > 
0(1 TeV), the bound on the lifetime becomes even more stringent (cf. [9J). We do not consider a late time 
entropy production in this paper, which is an another possible way to avoid these BBN constraints |10!. 
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strong constraints from hadronic showers [9]. Only a sneutrino NLSP could be marginally 
acceptable also with longer lifetimes, and therefore larger gravitino mass, in the region 
where the hadronic branching ratio of the decay is below 10~^ 

On the other hand, standard thermal leptogenesis [12], an attractive model for baryo- 
genesis, needs a large reheating temperature in the early universe, Tr > 10^ GeV 
(cf. [TSfT^ ). This reheating temperature implies 77^3/2 ^ 5 GeV for a gluino mass of 
rrig = 500 GeV in order to avoid overclosure of the universe due to thermal gravitino 
producfou [imsi The lower bouud on the grav,ti„o mass scales as mj;? ~ 

All these cosmological problems are automatically solved without any fine tuning 
of parameters in the case of a small breaking of R-parity, as given in Eq. ([2]), with a 
gravitino LSP. The NLSP lifetime becomes sufficiently short for A, A' > 10~^^, 



Therefore, primordial nucleosynthesis, thermal leptogenesis and gravitino dark matter 
are naturally consistent for 10~^^ < A, A' < 10^'' and 171^/2 ^ 5 GeV. This is the main 
point of this paper. 

The paper is organised as follows. In Sec. 2 we present a model where R-parity 
breaking is tied to B-L breaking, yielding the needed small R-parity breaking couplings. 
Sec. 3 deals with constraints from neutrino masses. Sec. 4 deals with implications for 
cosmology and collider physics. The results are discussed in Sec. 5. 

2 R-Parity Breaking and B-L Breaking 
2.1 A Model of R-Parity Breaking 

We consider a supersymmetric extension of the standard model whose symmetry group 
G includes U{1)b-l and R-invariance, 



Three quark-lepton generations can be grouped into the SU{5) representations 10, = 
(g, u'^, e^)i, 5j = {df^, l)i and 1 = vl, which together form 16-plets of 5*0(10). In addition, 
we have two Higgs doublets Hu and Hd-, two standard model singlets A^^ and A^, and three 

^We use the perturbative result for the gravitino production rate to leading order in the strong gauge 
coupling g. Since g and also the thermal gluon mass are large, the perturbative expansion is problematic 
[T5] . The uncertainty due to higher orders in g and nonpcrturbative effects is 0(1). Possible effects due 
to thermal masses are also 0(1) ,17|. 




G = SU{3) X SU{2) X f/(l)y X U{1)b-l x U{1) 



(5) 
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SO{10) singlets X, $ and Z. The two Higgs doublets are contained in 5- and 5-plets of 
SU{5), which we shall also denote as and H^, respectively. N'^ and are contained 
in 16 and 16 of 5*0(10), which fixes their B-L charge to be +1 and —1, respectively. X, 
$ and Z have B-L charge zero. This set of fields is familiar from 5*0(10) orbifold GUTs 
(cf. [H]): matter fields form complete 50(10) representations, whereas fields which break 
SU{2) X U{1)y and U{1)b-l appear as 'split' multiplets. For simplicity, we shall use in 
the following often SU (5) notation also for the Higgs multiplets. 

The matter sector of the superpotential has the usual form 



Wm = /iSf 10aO,//„ + hf^d^jHa + h^:;%l,H, + -L/i(;)l4^.iV2 , (6) 



where Mp = 2.4 x 10^® GeV is the Planck mass. The expectation values of the Higgs 
multiplets Hu and Hd generate Dirac masses of quarks and leptons, whereas the expec- 
tation value of the singlet Higgs field N generates the Majorana mass matrix of the 
right-handed neutrinos 1,. The superpotential responsible for B-L breaking is chosen as 

Wb-l = X{NN'^ - $2) , (7) 

where unknown Yukawa couplings have been set equal to one. $ plays the role of a 
spectator field, which will finally be replaced by its expectation value, {^) = vb-l- 
Similarly, Z is a spectator field [f, which breaks supersymmetry and U{l)ji, (Z) = FzOO. 
The superpotential in Eqs. ([6]) and ([7]) is the most general one consistent with the 
R-charges listed in Table 1, up to higher order terms which we will discuss later. Note 
that the choice of a negative R-charge for N"^ forbids the dangerous superpotential terms 

liHdN' , B.S.-lOfciV^ , (8) 

which would yield too large bilinear mixings and too rapid proton decay, respectively. 
The expectation value of $ leads to the breaking of B — L, 

(N) = {N') = ($) = VB-L , (9) 

where the first equality is a consequence of the U{1)b-l D-term. This generates a Ma- 
jorana mass matrix M for the right-handed neutrinos with three large eigenvalues, with 
Ml < Ms < Mg. If the largest eigenvalue of /i^") is C(l), one has Mg ^ t;|_^/Mp. 
The heavy Majorana neutrinos can be integrated out yielding for the matter part of the 
superpotential 

Wm = hff lOdO.H^ + la^.Hd - i(/i(^)^/i(^)^),,(5,if„)(5,iJ„) , (10) 
with the familiar dimension-5 seesaw operator for light neutrino masses. 



■^For simplicity, we use a spectator chiral superfield to describe supersymmetry breaking. The field Z 
is not essential for the connection between R-parity breaking and B-L breaking discussed in this section. 
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Table 1: R-charges of matter fields, Higgs fields and 5*0(10) singlets. 



Since the field $ carries R-cliarge —1, the VEV ($) breaks R-parity, which is con- 
served by the VEV (Z). Thus, the breaking of B-L is tied to the breaking of R-parity. 
This is the key feature of the mechanism for R-parity breaking presented in this pa- 
pei0. The breaking of R-parity is transmitted to the low-energy degrees of freedom via 
higher-dimensional operators in the superpotential and the Kahler potential. The leading 
correction to the Kahler potential is 

5Ki = -r^iaiZ^ + a[Z)<^^N%Hu + -^{ciZ^ + c[Z)^N%Hu + h.c. . (11) 
IVip IVip 

Replacing the spectator fields Z and $, as well as A^'^ by their expectation values, one 
obtains the correction to the superpotential 

5Wi = piiQliHu , (12) 

with 

^. = o(»3/.), e = |;f-f . (13) 

where 7/13/2 = Fz/{\/^Mp) is the gravitino mass. Note that O can be increased or 
decreased by an appropriate choice of Yukawa couplings in Eqs. (E]) and (171). Eq. ( |T2l) is 
the familiar bilinear R-parity breaking term [3]. The correction to the Kahler potential 

5K, = ^Z^H.H^ + h.c. (14) 



yields the corresponding R-parity conserving term [2T] 

5Wo = f^H^Hu , /i = C(m3/2) . (15) 

Note that fi and fii are generated by operators of different mass dimension. Hence, their 
values may easily differ by one or two orders of magnitude, allowing for m3/2 and 

a gravitino LSP. 

To analyse the complete superpotential including the R-symmetry breaking terms, 
it is convenient to perform a rotation of the Higgs and lepton superfields, 

Ha = H',- ed[ , = + e,H'^ , (16) 



*For a recent discussion of the connection between B-L breaking and R-parity breaking in the context 
of string compactifications, see [19|20j . 
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where = fiiO/fi. In terms of the new fields the superpotential reads 



W = Wm + SWo + 6Wi 

= f,H',H^ + hffq,u^^H^ + hfdtq,H', + h^h[e'jH', (17) 

The mixing of Higgs and lepton superfields has induced trihnear R-parity breaking terms 
(9(e). As we will discuss in Sec. 3, the mixing terms induce vacuum expectation values 
for the scalar neutrinos that in turn induce mixing terms 0{e) of neutrinos with the 
neutralinos, and neutrino masses suppressed by 0{e^). 

It is remarkable that the potentially dangerous operator leading to proton decay 
is strongly suppressed compared to the trilinear terms 0(e) in Eq. (fTTl) . The leading 
operator is 

5W2 = -^u^d'd^N^<l>^X . (18) 

For global supersymmetry one has (X) = 0, which in supergravity is modified to (X) = 
0(777.3/2) §• One then obtains 

SW, oc ^^^^I^^u'^d^d'^ + . . . . (19) 

For A, A' satisfying Eq. ([21), the coefficient of the dangerous dimension-4 Ai? = 1 operator 
is much smaller than the upper bound from the proton lifetime [1] . 



2.2 Scale of B-L breaking and Thermal Leptogenesis 

The phenomenological viability of the model depends on the size of R-parity breaking 
mixings and therefore on the scale vb-l of R-parity breaking. An important constraint 
comes from baryogenesis. As already discussed in the introduction, the potential washout 
of a baryon asymmetry before the electroweak phase transition is avoided if the R-parity 
violating Yukawa couplings satisfy Ajj^, A-^^ < 10~^, which in turn implies: 

This is a sufficient condition, which can be relaxed for some fiavour structures [5]. 

As an illustration for possible scales of B-L breaking we use a model [22] for quark 
and lepton mass hierarchies based on a Froggatt-Nielsen fiavour symmetry. The 

^The VEV (X) — O {1712/2) also causes an additional contribution to the bilinear term via SW = 
(l/Mp'^)X<f>Af^5ji7„, which is comparable to those from SKi. 
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Table 2: Chiral charges: a = or 1, ano? Q <h < c < d. 



mass hierarchy is generated by the expectation value of a singlet field with charge 
Q(f, = —1 via nonrenormalizable interactions with a scale A = {(f))/T] > Agut, V — 0.06. 
Yukawa couplings and bilinear terms for SU (5) multiplets ipi with charge Qi scale like 

hij oc 1]^^+^^ , oc . (21) 

Charges Qi describing qualitatively the observed quark and lepton masses and mixings 
are listed in Table 2. The model also predicts the observed baryon asymmetry via lepto- 
genesis for the cases where a + d = 2. There are two, at low energies indistinguishable, 
consistent scales of B-L breaking: M3 ~ 10^^ GeV {a = b = 0, c = 1, d = 2) and 
M3 ~ 10^2 GeV (6 = c = 0, a = d = 1). For fii/jj = 1.0 . . . 0.01 these two cases lead to 
the R-parity breaking mixing parameters (cf. Eq. ( fT3l) ) 

(I) 4^ = 10-3 . . . 10-5 , (II) 4t = 10-*^ . . . 10-^ . (22) 

In the extreme case M3 ~ M2 ~ Mi ~ 10^° GeV without Froggatt-Nielsen symmetry, 
where leptogenesis may still work for an appropriate enhancement of the CP asymmetry, 
one has 

(III) Q = 10-^ . . 10-1° . (23) 

In the fiavour models (I) and (II) the RPV mixings are suppressed by ri^\ As we shall 
see in the following section, model (I) is inconsistent with the constraints from neutrino 
masses and baryogenesis washout; the models (II) and (III) are consistent with both 
constraints. 

The expected mass scale of right-handed neutrinos depends on the mechanism which 
breaks B-L. The expectation value of a field with lepton number L = 2 can generate 
heavy Majorana masses via renormalizable Yukawa couplings. With B-L broken at the 
GUT scale, and for Yukawa coupling 0{1) for the third family, one then obtains the 
canonical result M3 ~ vb~l ~ 10^^ GeV. On the other hand, if right-handed neutrino 
masses are generated via a nonrenormalizable dimension-5 operator and the expectation 
value of a field with L = 1, as in Eq. ([6]), one has instead M3 ~ f^_^/Mp ~ 10^^ GeV. 
This illustrates how the two mass scales for M3, which correspond to the two cases (I) 
and (II), respectively, might be obtained. 
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3 Neutrino Masses 



The model we are considering generates after supersymmetry breaking bilinear R-parity 
violating terms, Eq. (IT^ . and tiny R-parity violating Yukawa couplings, Eq. (HM . that 
we neglect in what follows. Scenarios with just bilinear R-parity violation have been 
thoroughly studied in the literature [23]. Here, we will limit ourselves to estimate the 
size of neutrino masses, following closely |24j . 

At the high-energy scale, the soft SUSY breaking Lagrangian reads 

-Csoft = mjjjHd\^+mljHX+ml\h\^ + {BHdH^ + BjiH^+mlHjiH*a + h.c.)^ (24) 

For the computation of neutrino masses we find convenient to work in the basis where 
the R-parity violating bilinear couplings in the superpotential are rotated away, = 0, 
through the field redefinition Eq. ( fT6l) . This choice of basis has the advantage that once 
the basis has been fixed at the high energy scale, the condition /x, = holds at any 
scale, and it is not necessary to redefine the basis again at low energies. We also choose 
the phases of the lepton doublets such that the are real. In this basis the soft SUSY 
breaking Lagrangian is given by 

-C,,ft = mj,,\H'f+mljHX+ml\l[\'HB'H',H^+B'J[^^^ (25) 

where 

= - eiRe(m2^J + 0{\e,\^) , 
B' = B + BiEi , 
B[ = Bi — Bei , 

^Ih', = ^Ih, + e.(m2 - ml^) + 0(|e,n . (26) 

Minimisation of the scalar potential yields non-vanishing vacuum expectation values for 
the neutral components of the Higgs doublets, as well as for the sneutrinos, 

Bj tan B + ml „, 3 

These vacuum expectation values induce mixings between neutrinos and gauginos, giving 
rise to one non-vanishing neutrino mass through an "electroweak" seesaw mechanism. 
The resulting neutrino mass is 

= \9lvl E ^ , (28) 
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where gz is the Z-boson gauge couphng, m^o are the neutrahno masses and Cza is the 
Zino component of the neutrahno Xa- In the foUowing estimates we wiU replace the sum 
over inverse neutrahno masses by the inverse of the characteristic SUSY breaking scale, 
1/m. 

The size of the neutrino masses depends crucially on the mechanism for supersym- 
metry breaking. Generically, one expects 

B[ ~ /ijm ~ ejm^, rnfirri = eim^ . (29) 

i d 

Then, using Eqs. f l27|) and fl28|) we find for ei, 62 < €3, 

(30) 

We can now insert the values of €3 for the different fiavour models considered in the 
previous section. Clearly, model (I), which has a high scale of B-L breaking, is excluded. 
In models (II) and (III) one has 63 < 10^^ and < 10~^, respectively. Here the neutrino 
mass terms induced by R-parity breaking are negligible and the baryogenesis constraint, 
Eq. fl20l) . is fulfilled for all values of tan/3. 



4 Cosmology and Collider Physics 
4.1 Gravitino Decay 

Since R-parity is broken, the gravitino is no more stable in our setting, but it still has 
a lifetime much longer than the age of the universe, since it is suppressed both by the 
Planck mass and the small R-parity breaking parameters. 

The two-body decay is determined by the mixing of the neutralinos with the neutri- 
nos. Neglecting the small neutrino masses, one has 

r«'3/. - 7^) = . (31) 

The photino-neutrino mixing can be approximated by (cf. Eq. fl27j) ) 
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U^i> I — ()z 



for €1,62 < 63, where we made the rough estimate 0.1/m for the weighted sum of neu- 
trahno masses and the coupling, taking into account that not all mixings can be maximal. 
Using Mp = 2.4 x 10^*^ GeV, one obtains for the gravitino lifetime [6] 

^-.4.10=.(^)-^(^)^(^)-\ (33) 
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The three-body decay is usually subdominant due to the phase-space and interme- 
diate heavy particle suppression. For the decay with intermediate heavy tr, neglecting 
all external masses in the phase space factor, we find 

r(z^3/2 - rnkh) = ( ^] , (34) 

where 

F{a) = / dx—^ ^ -— . 35 

The full expression has been obtained in [25] . 

In the case where only bilinear R-parity breaking is present, the Xijk couplings are 
generated from the Yukawa couplings as 

Kjk = (^ihfk ■ (36) 
Then the inverse partial width for the three-body decay, 

r„.3,..«,-'.2x:o3'.(^)-(^^)-^(^)\^)-' , ,3T) 

is much larger than the lifetime determined from the two-body decay, Eq. (!33l) . as long 
as the mixing between photino and neutrino is not unnaturally suppressed. 



4.2 Extragalactic Diffuse Gamma-Ray Emission 

A stringent astrophysical constraint for decaying gravitino dark matter is the mea- 
sured gamma-ray fiux |6]. Assuming that the gravitino constitutes the dominant 
component of dark matter, its decay into neutrino and photon gives rise to an 
extragalactic diffuse gamma-ray flux with a characteristic energy spectrum, cor- 
responding to a red shifted monochromatic line. A photon with measured en- 
ergy £' = 7713/2/(2(1-1-2;)) has been emitted at the comoving distance x{z)y with 



dx/dz = (1 -|- 2;)^^/^/(ao-ffoY ^Af(l + '^(l + z)^^)). Here and Hq are the present scale 
factor and Hubble parameter, respectively, and k = Q/^/Qm — 3, with Q\ + Qm = 1, 
assuming a flat universe. For the photon flux one obtains, for T3/2 ^ H^^, 

E^^ = c,|i + .(£^) 1 (£^1 o[i-i^], {31 



dE 




with 



= . ""T^u. = 10-^ (cm^str s)-GeV f^)" ; (39) 



" 8nr,/,Hon\{' VlO-s 
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here r3/2 is given by Eq. (l33l) . and we have taken the gravitino density equal to the Cold 
Dark Matter density as ^3/2/^^ = 0.1, = 1.05 x IQ-^GeVcm-^, = 0.25 and 
Ho = h 100 km s"^ Mpc"^ with h = 0.73 [26]. 

In addition to the extragalactic signal one also expects a sharp line from the halo 
of our galaxy with an intensity comparable to the extragalactic signal and strong 
anisotropy [27] . We have in fact from the decay of halo gravitinos 



dE ^ \ 1713/2) 



E^=^=D^5(l-f^) , (40) 

where 



D, = C;-i^^l Phaio{l)dl- (41) 



The ratio D^/C^ is given only by cosmological constants and the halo dark matter density 
integrated along the line of sight, so the intensity and angular distribution of the halo 
signal is very sensitive to the distribution of the dark matter in the Milky Way. It is 
surprising that for typical halo models, such number is of order unity [23 cind shows 
moderate angular dependence if one excludes the galactic centre region. The anisotropic 
part of the halo signal may be partially hidden in the diffuse galactic 7-ray emission due 
to conventional astrophysical processes. We expect therefore the isotropic signal in the 
extragalactic 7-ray flux to be a combination of both the continuum spectrum in Eq. (l38l) 
and part of the halo line in Eq. ( HO|l . 

Assuming that one understands the diffuse galactic 7-ray flux, one can extract from 
the EGRET data the extragalactic diffuse component. The first analysis of Sreekumar 
et al. [28] gave an extragalactic flux described by the power law 

E2^ = 1.37x lO-M^— ) (cm^str s)-iGeV (42) 
dE VI GeV/ 

in the energy range 50 MeV-10 GeV. A non-observation of a 7-ray line can then be used 
to constrain the allowed gravitino mass and lifetime [6J. Assuming the gravitinos to make 
up all the Cold Dark Matter density, and taking a 3a upper bound on the flux above 100 
MeV corresponding to 2.23 x 10~^(cm^str s)~^GeV p8], we can have directly a rough 
bound on the gravitino lifetime from as 

r3/2 > 4 X lO^^s . (43) 

The more recent analysis of the EGRET data [29] shows in the 50 MeV - 2 GeV range 
a power law behaviour, but a clear excess between 2 GeV and 10 GeV. The maximal 
flux allowed by the data taking into account the model dependence and systematic 
errors is not very far from the one obtained in the old analysis, in fact the integrated 
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flux between 0.1-10 GeV is given as (11.1 ± 0.1) x 10 ^cm ^ str ^ s ^ compared to 
(14.5 ±0.5) X 10-6cm-2 str^^ s^^ [29]. 

This is precisely the energy range where, based on our lower bound on the gravitino 
mass of 5 GeV, one may expect a gravitino signal. It is very remarkable that also the 
measured flux corresponds to the expectation of the model for R-parity and B-L breaking 
discussed in Sec. 2 as can be seen from the bound Eq. (H3|) . On the other hand we would 
expect also an anisotropic flux from the halo component that EGRET does not resolve 
probably due to the galactic background, which is an order of magnitude larger than the 
extracted extragalactic signal. 

The excess in the extragalactic 7-ray flux above 2 GeV from the EGRET data [29] 
has also been related to the annihilation of heavy neutralinos in the galactic halo [30] . 
Due to the current limitations in the determination of the diffuse galactic 7-ray emission 
[31] theoretical interpretations of the EGRET excess remain uncertain at present. Clarifl- 
cation can be expected from the Gamma Ray Large Area Space Telescope (GLAST) [32] . 
to be launched this fall, that aims to improve by a factor 30-50 the sensitivity of the 
EGRET satellite to the diffuse gamma ray flux in the range 20 MeV-10 GeV. 

Another constraint comes from the neutrino flux. In the energy range of interest, 
from about 1 GeV to 1 TeV, the extraterrestrial neutrino flux is constrained by the 
flux of upward- going muons measured by the 1MB experiment, that does not show any 
discrepancy with respect to the expected neutrino flux from cosmic ray interactions 
in the atmosphere. The requirement that the neutrino flux from gravitino decay does 
not exceed the observed flux, translates into a lower bound on the gravitino lifetime, 
which is roughly r3/2 ^ 6 x 10^^ s for m^i2 = 1 TeV, and becomes weaker for smaller 
gravitino masses [33j. This bound is clearly consistent with a signal in the EGRET data, 
as discussed above. 

4.3 Collider Signatures 

The collider signatures depend on the nature of the NLSP. Here we consider the cases 
that the NLSP is the lightest stau or the lightest neutralino. 

The lightest stau, that we assume mainly right-handed, decays through tr 
T z/^,yU Ur- Ou the other hand, the small left-handed component of the stau mass eigen- 
state can trigger a decay into two jets through fi hH, provided the process is kine- 
matically open. The hadronic decays are enhanced compared to the leptonic decays by 
the larger bottom Yukawa coupling and by the colour factor, but are usually suppressed 
by the small left-right mixing. 

If the leptonic decay channel is the dominant mode, the decay length can be approx- 
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imated by 

cri- ~ 30 cm ( " ( ^) " f ^ V\ (44) 
V200GeVy VlO-V V 10 / 

It is intriguing that the sufficient condition to avoid the erasure of the baryon asymmetry, 

Eq. fl20|) . imphes the observation of a displaced stau vertex at future colhders, more than 

3mm away from the beam axis for 62 < 10~^. In the particular case of the flavour model 

(II) discussed in Sec. 2.2, 62 ~ 6 x 10^^, one has a spectacular signal consisting on a 

heavily ionising charged track of length ~ 0.8 m, followed by a muon track or a jet and 

missing energy, corresponding to r ^ /xz/,- or f — > rz/^, respectively. 

If the hadronic channel tl hH is the dominant mode, the decay length is given by 

/ nif { £3 ( tan/5\ ^ / cosOr 



-""-^■^-Isod^j lio^j [-ifj lirfj • («) 

where 6r denotes the mixing angle of the staus. This channel also yields a very unique 
signature at colliders, consisting of a heavily ionising charged track followed by two jets. 

These characteristic signatures would allow to distinguish at colliders our scenario 
from the case with conserved R-parity where the decay f — >• 1P3/2T leads to (cf. [34j) 

cr?/' ~ 40 cm ( ^] ' ( (46) 
UkeVy l200 GeVy " ^ ' 

Hence, for a gravitino mass 1713/2 ^ 0(10 keV), the decay length of the lightest stau 
is shorter than C(10 m), and would therefore decay inside the detector into tau and 
gravitino. The experimental signature for this process would be identical to the decay 
f Tv^. However, the scenario with R-parity violation also predicts the decay f /iz/^, 
with very similar branching ratio due to SU{2) invariance. Although this signature could 
be mimicked by a scenario with conserved R-parity if lepton flavour is violated, through 
the decay r — > /i'03/2, large branching ratios are precluded from present bounds on 
flavour violation [52]. In consequence, the observation of a comparable number of tau 
and muon events in stau decays would constitute a signature for the scenario with R- 
parity violation. Also, the observation of a stau decaying into two jets would undoubtedly 
point to the scenario with R-parity violation. 

On the other hand, if the lightest neutralino is the NLSP, it decays through Xi ~^ 
T^W^ [SB], or through Xi ~^ bb'^u if the former decay channel is kinematically 
closed. The corresponding decay lengths can be approximated by 

c.---' ~ 20 cm ( ) " (J^) -'(^J!^\\ (47) 

- . i^)' (^rc-^v . (48) 
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Again, this scenario can be easily discriminated at future colliders from the scenario 
with conserved R-parity. In this case, the neutralino decays into gravitino and photon 
[38] with decay length 

For a gravitino mass 7713/2 ^ (9(10 keV) the neutralino would decay inside the detector 
producing an energetic photon and missing energy, which is clearly distinguishable from 
the signals in the R-parity violating scenario that in general involve jets. 



4.4 Microscopic Determination of the Planck Mass 

Recently, a method has been proposed for the microscopic determination of the Planck 
mass in collider experiments [31], providing a direct test of supergravity. The method 
requires a very long lived stau NLSP which decays mostly into tau and gravitino, which 
is difficult to reconcile with recent constraints from BBN [7,.9|, unless there is a late-time 
entropy production [Tn|. In the picture proposed in this letter, where primordial nucle- 
osynthesis, thermal leptogenesis and dark matter are naturally consistent, this method 
cannot be pursued, as the stau decays predominantly in the R-parity violating channel 
into charged lepton and neutrino. 

Nevertheless, from a gravitino signal in the diffuse 7-ray flux and the width for the 
stau decay into two jets, one can still obtain a microscopic estimate of the Planck mass. 
The gravitino mass is given by the maximal energy of the photon, = 2E^, and the 
gravitino lifetime can be determined from the photon flux, Eqs. ( !38p0!) . Then, using the 
expression for the gravitino decay rate, Eq. (13T!) . one can rewrite the Planck mass in 
terms of the gravitino mass, lifetime and photino-neutrino mixing as 

= 2.6 X 10'« GeV ( ) I— ^) ISfl (50) 




where |f/^iy| is related to the decay rate of the stau into two jetq^l- We can cast the 
dependence on the decay rate as a dependence on the decay length of the stau in this 
channel, yielding 

lU- P ~ 10-- ('-^V ( ^ V C-^X ('-^X (51) 
l^^'^' 1.4m UoOGeV 10 0.1 ' 



^Note that the decay rate of the stau into leptons depends on 62 whereas \U^iy\ depends on the 
sneutrino VEV and therefore mainly on £3 for the hierarchical case. 
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The measurement of the decay length of the stau in the hadronic channel, com- 
plemented with additional information about supersymmetry breaking parameters, can 
provide a determination of \U^u\ - The measurement of the photon energy and the photon 
flux in the diffuse 7-ray background then gives the gravitino mass and lifetime and, using 
Eq. fl50l) . an estimate of the Planck mass. 

5 Conclusions 

On theoretical grounds, theories with and without R-parity are on equal footing. In this 
paper we have presented a simple model where R-parity is not conserved and its violation 
is connected to the scale of B-L breaking. One can then have R-parity violating couplings 
that are small enough to be consistent with baryogenesis and gravitino dark matter, yet 
large enough to allow for the NLSP decay before nucleosynthesis. For gravitino masses 
above 5 GeV one obtains a cosmological history consistent with thermal leptogenesis, 
thermally produced gravitino dark matter and primordial nucleosynthesis. 

Relic gravitino decays into neutrino and photon yield a diffuse halo and extragalactic 
7-ray flux which depends on the R-parity violating Yukawa couplings. It is remarkable 
that for a gravitino mass ^3/2 = C(10) GeV, the predicted photon flux could be part 
of the apparent excess in the extragalactic diffuse 7-ray flux obtained from the EGRET 
data. However, given the current uncertainties in the determination of the diffuse galactic 
7-ray emission, this consistency may be accidental. Unequivocal evidence for decaying 
gravitino dark matter could come from the results of GLAST. 

The flavour dependent pattern of R-parity breaking can give striking signatures at 
the LHC, in particular a vertex of the NLSP, that is significantly displaced from the 
beam axis. Together with the measurement of supersymmetry breaking parameters at 
the LHC, the observation of a redshifted photon spectral line from gravitino decay by 
GLAST can allow a microscopic determination of the Planck mass. In the less optimistic 
case where the R-parity breaking Yukawa couplings are near to their lower bound, astro- 
physical detection will be very challenging whereas signals hinting at R-parity breaking 
and gravitino dark matter could still come from stau decays, as in the case of R-parity 
conservation. 
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